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It is shown that the critical temperature for spin-triplet, p-wave superconductivity mediated by spin 
fluctuations is generically much higher in a Heisenberg ferromagnetic phase than in a paramagnetic 
one, due to the couphng of the magnons to the longitudinal magnetic susceptibility. Together with 
the tendency of the low-temperature ferromagnetic transition in very clean Heisenberg magnets to 
be of first order, this qualitatively explains the phase diagram recently observed in UGe2. 



It has long been known that, in principle, the exchange 
of magnetic fluctuations between electrons can induce su- 
perconductivity [1]. Magnetic fluctuations become large 
in the vicinity of continuous magnetic phase transitions, 
which makes nearly ferromagnetic materials, or ferro- 
magnets with a low Curie temperature, natural candi- 
dates for this phenomenon. In contrast to the much more 
common phonon-exchange case, which usually leads to 
electron pairing of spin-singlet, s-wave nature, the mag- 
netically mediated pairing is strongest in the spin-triplet, 
p-wave channel, p-wave superconductivity is very sensi- 
tive to nonmagnetic impurities and therefore can be ex- 
pected only in extremely pure samples. The combined re- 
quirements of high purity, low temperatures, and vicinity 
to a ferromagnetic transition severely restrict the number 
of promising materials. Indeed, until recently there were 
no convincing simple examples of magnetically induced 
superconductivity, and the paramagnon interpretation of 
superfluid ■^He [1,2] was considered the best example of 
pairing by exchange of magnetic fluctuations. 

This situation has recently changed, due to the obser- 
vation of the coexistence of ferromagnetism and super- 
conductivity in UGe2 [3]. In contrast to other uranium 
compounds, UGe2 has more in common with classic d- 
electron ferromagnets, like Fe, Co, and Ni, than with 
heavy-fermion systems. The persistence of ferromagnetic 
order within the superconducting phase has been ascer- 
tained by means of neutron scattering, and the itinerant 
ferromagnetism and the superconductivity are believed 
to arise from the same electrons [3]. Since superconduc- 
tivity in the presence of ferromagnetism must be of spin- 
triplet type, magnetically induced pairing is an obvious 
candidate for the observed superconductivity, although a 
phonon mechanism has also been proposed [4] . 

The nature of the phase diagram reported in Ref. [3] 
is, however, not obviously consistent with existing models 
of spin fluctuation induced superconductivity, see Fig. 1. 
Fay and Appel [5] have calculated the superconducting 
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FIG. 1. Schematic phase diagram showing the paramag- 
netic (PM), ferromagnetic (FM), and superconducting phases 
(SC) in a temperature (T) - control parameter (CP) plane, 
(a) shows the qualitative prediction of paramagnon theory [5], 
and (b) qualitatively shows the phase diagram as observed in 
UGe2 [3] and explained by the theory presented here. In Ref. 
[3] , hydrostatic pressure serves as CP. 



Tc for a p-wave, equal-spin pairing state in both the para- 
magnetic (PM) and ferromagnetic (FM) phases close to a 
continuous magnetic transition. Using a McMillan-type 
formula, they found values of Tc on either side of the tran- 
sition that are within 20% of one another. Their shape of 
Tc as a function of the distance t from the magnetic tran- 
sition is very similar to that obtained by Levin and Vails 
[2], who solved the Eliashberg equations numerically in 
the PM phase, although the absolute values of Tc are 
smaller in the McMillan approximation. More recently, 
Roussev and Millis [6] have obtained similar results in 
the PM phase. Contrary to this theoretical expectation 
of a superconducting phase diagram that is roughly sym- 
metrical with respect to the magnetic phase boundary. 
Fig. 1(a), the authors of Ref. [3] observed superconduc- 
tivity at temperatures up to about 500 mK within the 
FM phase only. Fig. 1(b). Qualitatively the same phase 
diagram has very recently been observed in ZrZn2 [7]. 
Since the spin fluctuations become large on cither side of 
the magnetic transition, it seems hard to reconcile this 
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experimental result with paramagnon theory [8]. 

In this Letter we show that the observed phase dia- 
gram can nevertheless be understood in these terms. The 
key lies in the existence of spin waves or magnons in the 
FM phase, which couple to the longitudinal susceptibil- 
ity and contribute a mode-mode coupling term to the 
latter that has no analog in the PM phase. We will see 
that this produces a superconducting transition temper- 
ature which under reasonable assumptions can easily be 
50 times larger than in the PM phase. 

We have included this effect in a McMillan-type Tc 
calculation similar to the one in Ref. [5] A representa- 
tive result of our analysis is shown in Fig. 2. The solid 
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FIG. 2. Superconducting Tc (solid curve, left scale) as a 
function of the distance from the critical point t, and the mag- 
netization A4. The dashed line (right seale) shows Tc in the 
PM phase scaled by a factor of 50, and the dotted curve (right 
scale) is the result in the FM phase without the mode-mode 
coupling effect. See the text for further explanation. 

curve (left-hand scale) represents the superconducting Tc 
as a function of the dimensionless distance t from the 
FM critical point. Also shown is the magnetization M 
in the FM phase, in units of the saturation magnetiza- 
tion fJ^Bn, with n the electron rmmbcr density and fxs 
the Bohr magneton. is measured in units of a charac- 
teristic temperature To that is given by either the Fermi 
temperature or a band width, depending on the model 
considered. The dashed curve shows the result in the PM 
phase scaled by a factor of 50 (right-hand scale), and the 
dotted curve in the FM phase (also scaled by a factor 
of 50, right-hand scale) represents the result that is ob- 
tained upon neglecting the mode-mode coupling effect. 
Notice that the maximum Tc in the FM phase is more 
than 50 times higher than in the PM phase. This relative 
difference between Tc in the two phases is the important 
result of our analysis. The absolute values should not be 
taken very seriously, as calculating Tc is notoriously diffi- 
cult and our simple mean-field treatment is certainly not 
adequate for this purpose. However, the relative compar- 
ison we expect to be reliable. It predicts a pronounced 
asymmetry between the PM and FM phases, which in 
the case of UGe2 means that superconductivity in the 



PM phase should not be expected at temperatures above 
at most 10 mK, in agreement with the experiment. Also 
of interest is the fact that at low temperatures the FM 
transition in very clean itinerant Heisenberg systems is 
generically of first order, as has been predicted theoret- 
ically [9] and is indeed observed in UGe2 [3] as well as 
in MnSi [10]. For the purpose of our discussion this sim- 
ply means that values oi \t\ smaller than some minimum 
value are not experimentally accessible. 

In the remainder of this Letter we sketch the theoret- 
ical analysis that has led to these results. For an order 
parameter (OP) field, we choose !F{x,y) = i()-\{x)'i()-\{y), 
with ipcr{x) a-n electronic field with spin index a and 
space-time index x [11]. The OP, i.e. the expectation 
value {J-{x, y)) = F{x — y), is the anomalous Green func- 
tion. The orbital symmetry of the OP is still unspecified, 
we will later choose the p-wave case. 

We have derived coupled equations of motion for F 
and the normal Green hmction, G, that lead to a loop 
expansion for the equation of state [12]. Our model is a 
microscopic action S with a free-electron part. So, and 
spin-singlet and spin-triplet interaction terms, 

SL = yJ dxl^sixf, S^^, = =^J dx[ncixf. (1) 

Here ns{x) and nc{x) are the electronic spin and charge 
density fields, respectively, and Ft and F^ are the spin- 
triplet and spin-singlet interaction amplitudes. We as- 
sume that screening has been built into the starting ac- 
tion, so the interaction amplitudes are simply numbers. 
By putting F^ = Fj one obtains the Hubbard model con- 
sidered in Ref. [5]. The magnetic equation of state we 
treat in zero-loop approximation. The superconducting 
equation of state needs to be calculated in one-loop ap- 
proximation in order to capture the spin-fluctuation in- 
duced pairing. It takes the form of linearized strong- 
coupling equations that are similar to those in Ref. 
[5]. These equations can be rewritten as an eigenvalue 
problem, which can then be solved numerically, using 
some theory for the (para)magnon propagators as input. 
This is the established procedure to calculate the criti- 
cal temperatin-e for phonon-mediated superconductivity 
[13], and it has been employed in the case of magnetically 
induced superconductivity or superfluidity in Refs. [2,6]. 

Even with a complete numerical solution of the strong- 
coupling equations, the superconducting Tc is notoriously 
hard to predict. This holds a fortiori in the case of mag- 
netically mediated superconductivity since (1) there is 
much less experimental information about the param- 
agnon propagator, that could be used as input, than 
about phonon spectra, and (2) there is no analog of 
Migdal's theorem. Our ambition here is therefore not 
to calculate but rather to perform a relative compar- 
ison of Tc values in the PM and FM phases, respectively. 
For this purpose, a simple McMillan-type approximation 
for Tc [13] suffices. We obtain 
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Te = To{t) exp [-(1 + 4 + 24)/df ; 



(2) 



Here To{t) is a temperature scale that will be specified 



below. Specializing to the p-wave case, the d^'^ read 
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kp{kp) are the Fermi wavenumbers for the up (down)- 
spin Fermi surface, and Np is the density of states at 
the up-spin Fermi surface. In the PM phase, kp — 
4 = kp- DL,T{q) are the longitudinal and transverse 
(para)magnon propagators. They are related to the elec- 
tronic spin susceptibility x via Dl,t{q) = Xl,t{q)/'2.Nf, 
with Np the density of states at the Fermi level in the 
PM phase. We use the expressions that were derived in 
Ref. [14], with one crucial modification that we will dis- 
cuss below. Prom that paper we obtain in the PM phase, 
and in the limit of small wavenumbers, 



£>i,T(g, iO) = l/[i + hL,T{q/2kFf 



(4) 



In the Gaussian approximation of Ref. [14], bL = br = 
1/3. However, there is no reason to prefer this Gaussian 
approximation over any other approximation scheme. 
The functional form of the long-wavelength expression, 
Eq. (4), on the other hand, is generic. We therefore adopt 
Eq. (4) with bL,T arbitrary coefficients of 0(1). By the 
same reasoning, we have in the FM phase, in the limit of 
long wavelengths and small frequencies. 



DL{q, iO) = l/[5|t|/4 + bLiq/^kpf] , 
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, (5b) 



with A the Stoner band splitting. For < A < nPt, A 
is related to the magnetization M by M = ^eA/Ft. 

Two comments: (1) In a strict long-wavelength expan- 
sion of the propagators from Ref. [14] the 6l and bx in 
Eqs. (5a,5b) become magnetization dependent. We ig- 
nore this efi'ect and use the same values as in the PM 
phase. We have compared this approximation against 
using the full propagators from Ref. [14], see below. (2) 
The factor of 5/4 in Eq. (5a) arises since we keep the par- 
ticle number density fixed, as is the case experimentally, 
rather than the chemical potential, sec Ref. [5] . 

We now consider the longitudinal magnetic propaga- 
tor in the FM phase in more detail. In a Heisenberg 



ferromagnct, the transverse spin waves or magnons cou- 
ple to the longitudinal susceptibility xl- This efi'ect 
is most easily demonstrated within a nonlinear sigma- 
modcl description of the ferromagnct, which treats the 
order parameter M as a vector of fixed length M, and 
parameterizes it as M = M {a{x),-K\{x),-K2{x)) with 
cr^ + TT^ + 7r| = 1, M the magnetization, and x a space- 
time index [15,16]. The diagonal part of the TTi-TTj prop- 
agator, (TTiTTi) = {M'^ /2Np) Dt, is proportional to the 
transverse propagator Dt, and the off-diagonal part has 
been calculated in Ref. [16]. The longitudinal propaga- 
tor, Dl = {M^ /2NF){a{x)a{y)), can be expanded in a 
series of 7r-correlation functions. 



{a{x)aiy)) = l-2{7r,Xa:)Tr\x)) 



(6) 



where repeated indices are summed over. At one-loop 
order, the term of order tt* yields the diagram shown 
in Fig. 3. Notice that the sigma model, which neglects 
all longitudinal fiuctuations, replaces the external legs by 
constants. Power counting shows that at nonzero tem- 




FIG. 3. Mode-mode coupling contribution to the longitu- 
dinal (L) propagator Dl from the transverse (T) ones. 

perature, and for dimensions d < 4, this contribution 
causes the homogeneous longitudinal susceptibility to di- 
verge everywhere in the FM phase [17]. More generally, 
this one-loop contribution, together with the zero-loop 
one, Eq. (5a), yields a hmctional form for Dj^ in the 
FM phase that is exact at small wavenumbers. This di- 
agram has no analog in the PM phase, while all other 
renormalizations of the propagators will give compara- 
ble contributions in both the PM and FM phases. It 
is therefore reasonable to calculate Tc based on a one- 
loop approximation in the FM phase, and compare it to 
a zero-loop calculation in the PM phase. We have used 
Eq. (5b) for the internal propagators in Fig. 3. Since the 
coupling constants involve a wavenumber integral, Eqs. 
(3), we also need to go beyond the sigma model and keep 
the wavenumber dependence of the external ones. For 
computational simplicity, we have modeled the external 
legs by replacing Eq. (5a) with a step function that cuts 
off the momentum integral at k/2kp — yjh\t\/ Ab^. With 
these approximations, the momentum integral in Fig. 3 
can be done analytically, leaving the frequency sum to 
be performed numerically. The result, and the resulting 
contribution to and , depend on the temperature, 
so Eq. (2) now needs to be solved self-consistently for Tc. 

We still need to specify the temperature scale TQ{t). 
Following Ref. [5], we use the prefactor of \t\ in Eqs. (4,5a) 
as a rough measure of the magnetic excitation energy. 
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To{t)=To[eit)t + e{-t)5\t\/4] , 



with To a microscopic temperature scale that is related 
to the Fermi temperature (for free electrons) or a band 
width (for band electrons). This qualitatively reflects 
the suppression of the superconducting Tc near the FM 
transition due to effective mass effects [2,5,6]. 

We are now in a position to choose parameters and cal- 
culate explicit results. We put Tg = Tt [5]; other reason- 
able choices yield similar results. Let us first ignore the 
mode- mode coupling contribution to and • We have 
performed the calculation both with the full propagators 
from Ref. [14] and with the schematic Landau propaga- 
tors, Eqs. (4) and (5a). With 6l = 0.23, br = 0.4 the 
two results arc within 10% of one another, and also very 
similar to those obtained by Fay and Appel [5] . We then 
use these values of bL,T to calculate the mode-mode cou- 
pling contribution, and solve the Tc-equation. The result 
is shown in Fig. 1 and has been discussed above. We have 
also explored the effect of varying the parameters bL,T- 
With 6l = &T = 1 we obtain the result shown in Fig. 
4. The (unphysical) zero-loop result in the FM phase is 
very sensitive to the parameters, while the enhancement 
of the (physical) one-loop result over the Tc in the PM 
phase is rather robust. However, the position of the max- 
imum of Tc changes compared to Fig. 1; it now occurs 
at the point where the magnetization reaches its satura- 
tion value. The reason is as follows. As one approaches 
the magnetization saturation point from low magnetiza- 
tion values, the transverse coupling constant c?g vanishes, 
and remains zero in the saturated region. Effectively, the 
Heisenberg system turns into the Ising model discussed 
in Ref. [6]. If the longitudinal coupling constant cif still 
has a substantial value at that point, then this leads to 
an increase in Tc- This is a very strong effect in the zero- 
loop contribution, see Fig. 4, and the effect qualitatively 
survives in the one-loop result. If, however, is already 
very small, then dp going to zero has only a small effect 
on Tc, as is the case in Fig. 1. Which of these two cases 
is realized depends on the parameter values. We finally 
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FIG. 4. Same as Fig. 2, but for different parameter values 
(see the text). 



(7) mention that the first order nature of the magnetic tran- 
sition [9,3] adds another mechanism for suppressing Tc in 
the PM phase: For a suflaciently strong first order tran- 
sition, and if the case shown in Fig. 4 is realized, then 
the effective t may be large enough everywhere for the 
system to miss the maximum of Tc in the PM phase, but 
not in the FM phase. 
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